Theory of mesoscopic magnetism in meta-materials 
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We provide a rigorous theoretical basis for the artificial magnetic activity of meta-materials near 
resonances. Our approach is a renormalization-based scheme that authorises a completely general 
theory. The major result is an explicit expression of the effective permeability, in terms of resonant 
frequencies. The theoretical results are checked numerically and we give applications of our theory 
to Left-Handed Media and to the solution of the Pokrovski-Efros paradox. 



PACS numbers: 73.20.Mf, 41.20.Jb, 42.70.Qs 

Photonic crystals are artificial devices, periodically 
structured, that exhibit photonic band gaps Di- 
electric photonic crystals are considered in the optical 
domain, but metallic ones (or wire mesh photonic crys- 
tals) are also studied in the microwave or TeraHz range 
y{. It has been well established Q that, below a 
cut frequency, wire mesh photonic crystals behave as 
if they were homogeneous with a negative, frequency- 
dependent, permittivity given by Scs = 1 — 2Trj/ [uo/cf , 
where 7 = log(d/r) (here d is the period of the crystal 
and r the radius of the wires) . For a fequency below 
LOp — -v/StF/c, the homogenized permittivity is negative 
and the propagation of waves is forbidden. The homog- 
enized permittivity represents the scattering behavior of 
the wire mesh photonic crystal for large enough wave- 
lengths, and explains why these structures present a pho- 
tonic band gap down to the null frequency (at least for 
infinitely conducting wires). Recently, Pendry and co- 
authors suggested that it was possible to design photonic 
crystals with non-magnetic materials that they possess 
an artificial magnetic activity and be described by 
an effective permeability. Basically, two geometries have 
been suggested: Split Ring Resonators and dielectric 
fibers with a large permittivity 0, • It is believed that 
with these geometries it is possible to obtain a negative 
permeability, and, by adding a wire mesh structure, to 
design a material with both a negative permittivity and a 
negative permeability. Materials with these characteris- 
tics do not seem to exist in nature, and therefore one tries 
to design them artificially (they are called "Left-Handed 
Materials"). They were studied theoretically long ago in 
a speculative and quite fascinating work @ by Veselago. 
He showed that they behave as if they had a negative 
index. Among other properties, Snell-Descartes law is 
reversed: at the interface between air and the material a 
beam is refracted on the same side of the normal. These 
ideas have motivated a lot of works, both experimentally 
and numerically (in particular in and also a lot of 



polemical issues [Toj . 

It seems however that, hitherto, there be no unified 
theoretical approach to this kind of effective behavior. 
The method generally followed consists in characteriz- 
ing the scattering matrix of a basic resonator by means 
of the £ and /i parameters, and then deriving the effec- 
tive parameters without taking into account the coupling 
between each resonator 11] . In the present work, we 
address this problem by using a renormalization group 
analysis, which gives us a deep insight into the phenom- 
ena and predicts that the existence of a negative /i is 
linked to internal resonances. In fact, the possibility of 
getting a negative permeability is very different from the 
possibility of getting a negative permittivity: while the 
negative e is obtained for low-frequencies (i.e. large wave- 
length with respect to the wires constituting the crystal) , 
the negative /i is obtained in the resonant domain, and 
only for a rather small interval of frequencies. In partic- 
ular, our approach explains the apparent paradox raised 
by Pokrovsky al. [l2|, that, by embedding wires in a 
medium with negative /x, one does not get a Left Handed 
Medium. It gives also a complete analysis of the effective 
properties of wire mesh photonic crystal, with a very high 
conductivity. Before going into the details of our study, 
we stress that we have tried here to make a bridge be- 
tween two domains that seem to be antagonist: that of 
the resonances and that of homogenization. 

In the following, we consider a 2D photonic crystal 
whose Wigner-Seitz cell Y is given in fig. ^ It is made of 
a dielectric rod (relative permittivity e;, cross-section D) 
embedded in a dielectric matrix Se- When the contrast 
between £i and Se is substantial enough, there appear Mie 
resonances into the hi ghe st index material. It was sug- 
gested some time ago |l3| that such internal resonances 
might result in the opening of forbidden gaps. Our point 
is to show that, near these resonances, the device behaves 
as if it had homogeneous electromagnetic parameters Sh 
and /i/j. Of course, for this situation to be physically 
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FIG. 1: The basic cell of the photonic crystal 

sounded, the resonant wavelengths should be larger than 
the period, otherwise the medium could not be described 
by homogeneous parameters. That is why we request, 
as in 0, that Si be much higher then Se- The method 
that we employ consists in changing ("renormalizing") 
the properties of the medium while keeping the relevant 
physical phenomena, i.e. the resonances, unchanged. To 
do so, we choose a small number r < 1, and we proceed 
to the following operation, denoted TZ: 

• We multiply the radius of the rods and the period 
by T, while maintening the domain O where the 
rods are contained constant (the number N of rods 
is increased as TV |f7|/T^). 

• We divide the permittivity e of the rods by (the 
optical diameter remains constant). 

The wave is p-polarized so the induction field reads 
B(x) = M(x)ez, but the vectorial form will prove use- 
ful for the analysis. We write B(x;7?.) and E(x;7^) the 
fields scattered by the renormalized structure. The point 
is to iterate this operation n times and study the limit of 
B(x;7?,"') and E(x;7?.") as n tends to infinity. In order 
to do so, we use a two-scale expansion of (E, B): 

B =Bo(x,x/T")-f r"Bi-t-... 

E = Eo(x,x/r") +t"Ei + ... ^' 

where the fields Eo,Bo depend on both the real space 
variable x (the global variable) and on the Wigner-Seitz 
cell variable y (the local variable). The fields are peri- 
odic with respect to y. Our point is to find the limit 
fields Eo,Bo. The local variable is in fact a hidden one: 
it is an internal degree of freedom. The true (observ- 
able) macroscopic fields (E^ , B/i) are the averages of the 
microscopic fields (Eo,Bo) over Y: 

B,,(x) = ^Bo(x,y)dy,E^(x) ==^Eo(x,y)dy. (2) 

Although we do not find it relevant to present all the 
mathematical details, we believe that it is important to 
offer the reader a general view of the method employed 
to get the limit fields. A complete and rigorous mathe- 
matical derivation can be found in 1151 



We analyze first the behavior of the fields with respect 
to the local variable. That is, we wish to describe the 
microscopic behavior of the fields with respect to their 
internal degrees of freedom. Using the expansion of 
the field, the V x • operator is transformed into: 

V X > Va; X • + T""Vy X • 

(we have to make explicit on what variables the deriva- 
tions operate because there are two sets of variables). 
Plugging these expressions into Maxwell system and 
identifying the terms that corresponds to the same power 
of r" we obtain the following system for the microscopic 
electric field: 

Vj, X Eo = on y , Vy • Eo = on F \ D (3) 

Besides: Eq = on D and Ei = on F \ D. This 
system is of electrostatic type: Eq does not depend on the 
microscopic induction field nor does it depend upon the 
wavelength. As a matter of fact, onY\D, Eg does not 
depend upon the variable y, as it can be deduced from 
system Q. Let us now turn to the magnetic field. The 
system satisfied by Bq is of an entirely different nature: 

Vy X Bo = -iujEiEi on Y , , 

Vy X El = icjBo on D ^ ' 

We have obtained a microscopic Maxwell system that 
describes the microscopic behavior of the fields. It can 
be seen that Ei gives indeed a first order expansion of the 
field inside D: it replaces Eq which is null there. Let us 
now use the fact that the fields are polarized. We write: 
Bo(x) = uo(x)ez. Plugging this expression into system 
Q shows that the magnetic field is independent of y on 
Y \ D. Next, by combining the equation is system Q), 
we find that: 

AyUQ + k^SiUQ = on 13 , uq = est onY \ D (5) 

We deduce from this system that the microscopic induc- 
tion field is linked to the macroscopic one by: uo(x, y) = 
(r7i(y)//^/i) u/i(x) where m satisfies: 

AyTO -I- k^e^m = on D , m = 1 on F \ D (6) 

and /ift, which shall be interpreted below as a relative 
permeability, is the mean value of to : — Jy "^(y)dy. 
Up to now, we have clarified what happens at the mi- 
croscopic scale. The point is now to derive the equations 
that are satisfied by the macroscopic fields. The propa- 
gation equations read, ior y E Y \ D: 

Vj; X Bq Vy X Bi = -iwEoeeEo 
Vj; X Eo + Vj, X El = iwBo 

In the first line, we recognize the Maxwell- Ampere equa- 
tion with the extra-term Vy x Bi. This term is homo- 
geneous to an electric displacement field, and it repre- 
sents the polarisation due to the presence of the scat- 
terers. Indeed, in the long wavelength regime, the emis- 
sion diagram of a fiber is that of a dipole (for the p- 
polarization) . As a consequence, the whole set of fibers 
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that constitutes the photonic crystal behaves as a set of 
coupled dipoles, producing a possibly anisotropic permit- 
tivity tensor. More precisely, as uq does not depend on 
y on F \ D, we obtain the following system satisfied by 
Mi: 

AyUi = onY\D , ^ -n-VyUo on D, (8) 

where n = (71,1,712) is the normal to D. This sys- 
tem implies a linear relation between ui and uq of the 
form: Vj^ui — 'P{y)VyUo where: 

(X I dwi dw2 \ 

and Wi satisfies: 

Aw, = on Y\D, = -n, on dD (10) 
on 

It is not diflicult to see [H [13 that Ah = Jy'P{y)dy 
is the inverse of the effective permittivity tensor eh{ = 
Aj^^) of the photonic crystal. The effective macroscopic 
equation can now be obtained by averaging system O 
on Y\D: 

^ ■{ej;^^{fij;^Uh)) + kW^O (11) 

The macroscopic behavior of the system is characterized 
by an effective permittivity tensor Eh and an effective 
permeability ■ This shows that the system exhibits 
an artificial magnetic activity. There are two huge differ- 
ences between the effective permittivity and the effective 
permeability: the permittivity can be a matrix, hence the 
medium can be anisotropic, whereas the effective perme- 
ability is always a scalar, therefore no anisotropic perme- 
ability can be obtained. Second, the permittivity is not 
frequency dependent, it is a static permittivity, whereas 
the permeability depends on the frequency. Let us give 
a closer look at the system of equations that defines the 
effective permeability ^h- As it stands in (|HJ it is just 
a partial differential equation problem. Under this form, 
its physical meaning does not appear clearly. To make it 
more explicit, let us recast it into an eigenvalue problem. 
This will help us understanding the underlying physics 
of what might look, at this stage, as a rather abstract re- 
sult. The system has a unique solution only if there is 
no function tfj such that 4' is null on Y\D and 4' satisfies 
the same Helmholtz equation on D. Otherwise m -\- tp 
would still be a solution of 10. Following spectral theory 
[T^ . we denote H = — £~^A and we look for functions $ 
satisfying the eigenvalue problem: 

^ ^0 onY\D , = k^'^ onD. (12) 

We get a set of eigenvalues k"^ and a set of eigenfunctions 
|$„). The physical meaning of these eigenvalues can be 
understood by going back to the unrenormalized initial 
fiber, with permittivity Si. This fiber alone possesses res- 
onant frequencies. They correspond to modes that are 



strongly localized inside the fiber. However, when there 
is a large number of fibers these resonances are slightly 
shifted due to coupling, and these resonances are further- 
more modified by the renormalization process. That is 
precisely what the eigenvalues of problem (|12|l are: the 
renormalized Mie frequencies of the fibre. 

For a given wavevector fc^ , we look for a solution m by 
expanding it on the basis |<&„), by noting that m — 1 is 
null on Y\D: m{y) = 1 + r7i„ . The coefficients 
are obtained by inserting this expansion in l^. We get, 
after averaging, the effective permeability — {l\m) 
under the form: 

We have obtained a completely general expression for the 
effective permeability. It relies on the cavity modes of 
the fibre only. In the vicinity of a resonance fc^, we have: 

^ ^ — kn |(*i'ri|l)r (k^ — k^iY^ which shows, in com- 
plete generality, that the permeability exhibits anoma- 
lous dispersion near the resonances, and becomes nega- 
tive there. It should also be noted that only the eigen- 
functions with non-zero mean value contribute. This is 
due to the fact that we have only kept the first order 
terms in the expansions JQ. 

Let us give an explicit computation in case of a square 
fiber. The derivation is rather straightforward, and 
follows closely that of the well-known TE modes in 
waveguides with square section. The eigenfunctions are 
^nm(y) — 2 sin{mryi) sin{rmTy2) and the corresponding 
eigenvalues are: fc^„ — ■n'^{n'^ -\- m?) The expansion of m 
on this basis leads to the following effective permeability: 

{n,m)odd ^ ^ 

where /c^„ = k'^^/a^Ei. Let us now turn to some nu- 
merical applications. First, we note that our analysis is 
supposed to work when there are Mie resonances at wave- 
lengths large with respect to the period of the crystal. 
This was the situation described in 7], where e = 200-|-5i. 
We choose this value for our own numerical computa- 
tions, the point being to test the vahdity of the theory. 
Using a rigorous diffraction code for gratings 0, we 
plot the transmission spectrum (dashed line in fig. ^ 
for a stack of 3 diffraction gratings made of square rods. 
There is a band gap for X/d between 8 and 12, due to a 
Mie resonance. In order to test our results, we plot the 
transmission spectrum of a homogeneous slab (solid line 
fig. 121) with parameters Eh = 1-7 (this value is obtained 
numerically from the resolution of problem (|10|l ) and /i^ 
given in (|14|) . We see in fig. |21that both curves fit very 
well, indicating that, although the wavelength is not that 
large, the whole photonic crystal behaves as a homoge- 
neous magnetic material. The discrepancy that is seen 
around X/d = 6.5 is due to the presence of a Mie reso- 
nance of null mean value that is not taken into account 
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FIG. 2: Real (bold lines) and imaginary (thin lines) parts of 
the transmission for the meta-material (solid lines) and the 
homogenized material (dashed line). 

in our theory. Only by expanding the fields to the second 
order could we incorporate this resonance in our global 
result. 

Finally, let us use these results to analyze recent prob- 
lem. In Pokrosky and al. showed that it was not 
possible to design a negative index medium by embedding 
metallic wires into a matrix with a negative fi, whereas 
the converse is possible. This can be explained in the 
following manner: the negative permittivity is obtained 
as a macroscopic effect, by which we mean that it is an 



interference effect and not an effect that takes place at 
the scale of the microscopic cell only. In much a different 
way, the negative permeability is obtained as a purely 
local effect, that happens at the scale of the microscopic 
cell. Therefore, for this effect to occur, no strong coupling 
between the fibers is requested, the coupling has just to 
be sufficient enough that the incident field can reach the 
fibers by tunnel effect. In our model, the propagation 
equation of the structure is obtained immediately by re- 
placing Ee by —Eg- Then, near the regions of anomalous 
dispersion, both parameters are negative and the prop- 
agation equation is the usual Helmholtz equation. Con- 
sequently, the field can propagate. On the contrary, for 
metallic wires in a medium with negative /x the propaga- 
tion equation reads T7,'T^: Au+k^ ii{1 — 2tt"/ / fik^)u — 0, 
which leads to evanescent waves. 

We have given in this work a theory of the mesocopic 
magnetism in metamaterials. We have shown that it was 
possible to give a homogenized description of a hetero- 
geneous device in the resonance domain. To do so, we 
have used a renormalization approach that shows that 
two scales should be distinguished: a microscopic one 
and a macroscopic one. We have shown that the artifi- 
cial, mesoscopic, magnetism is due to microscopic mag- 
netic moments induced by geometric resonances. So far, 
the analysis works for high permittivities, but we stress 
that there are inner resonances in gratings for much lower 
contrasts as well ^20]. Therefore, we do believe that the 
same physics can be found in the optical range of wave- 
lengths. 
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